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Bose-Einstein condensates confined in traps exhibit unique features which have been the object 
of extensive experimental and theoretical studies in the last few years. In this paper I will discuss 
some issues concerning the behaviour of the order parameter and the dynamic and superfluid effects 
exhibited by such systems. 



A. Introduction 



The aim of these lectures is to review some key features exhibited by Bose-Einstein condensation in trapped atomic 
gases, with special emphasis to the phenomena of superfluidity which have been the object of recent experimental 
and theoretical investigation and which are characterized in a non trivial way by the combined role of the phase of 
the order parameter and of two-body interactions. 
t-H , Let us start recalling that these systems are characterized, at sufficiently low temperatures, by a complex order 
^ ■ parameter which can be written as 

o\ ■ 

OV * (r,t) = AAuM)e iS ^ (1) 

where no =| | 2 is the so called condensate density and S is the phase. In general the condensate density does 
not coincide with the density n of the gas, because of the occurrence of quantum and thermal effects. However, in a 
very dilute gas, where the quantum depletion is negligible and for low temperatures, such that the thermal depletion 
is small, the condensate density coincides in practice with the total density. For this reason one often identifies 
n with n. For the same reason in these systems the condensate density can be measured with good precision 
through imaging techniques. Even at higher temperature the condensate density can be extracted experimentally 
with reasonable accuracy. In fact, due to the harmonic shape of the external confinement, the thermal component 
is spatially separated from the condensate so that the density profile is characterized by a typical bimodal structure 
allowing for systematic experimental analysis. The order parameter of trapped atomic gases then emerges as a crucial 
quantity not only from the conceptual point of view like in all many-body systems characterized by a broken symmetry, 
but also from the experimental one. This represents a key difference with respect to superfluid helium where the order 
parameter is a very hidden variable which can be measured only under very special conditions. The order parameter 
of a trapped gas can vary in space, time and temperature and characterizes in a direct way most of the physical 
!-h ' quantities that one measures in these systems. 

In the first part I will mainly discuss the properties of the modulus of the order parameter. From the theoretical 
point of view this quantity is interesting because even its shape at equilibrium is very sensitive to the presence of 
two-body interactions. In the second part I will give more emphasis on the role of the phase, a quantity which plays 
a crucial role in characterizing the coherence and superfluid features of our systems. Of course the modulus and the 
phase of the order parameter are deeply related quantities and in many cases their behaviour cannot be discussed in 
a separate way. Many of the results discussed in these notes are derived in more detail in ref. [lj to which we refer 
also for a more complete list of references. 

The discussion on Bose-Einstein condensation presented in these notes is stimulated by several features of funda- 
mental interest that are worth mentioning. In particular it is important to recall that 

- BEC gases in traps provide a unique opportunity to study the transition from the microscopic to the macroscpic 
world. 

- In these systems BEC shows up both in momentum and coordinate space. 

- BEC gases provide an almost ideal realization of classical matter waves, giving rise to new coherence phenomena. 

- Trapped gases are well suited to study rotational superfluid phenomena. 
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B. BEC is a phase transition 



In the presence of harmonic trapping 

Vext = \m (uj 2 x x 2 + uj 2 y y 2 + uj 2 z z 2 ) (2) 

the critical temperature for Bose-Einstein condensation of an ideal gas can be easily evaluated using the standard 
techniques of quantum statistical mechanics. To this purpose one has to define the thermodynamic limit. In the case 
of harmonic trapping the natural choice is to take N — > oo and 0Jh o — ► 0, with the combination Nu^ kept fixed. 
Here ujho = {^x^y^z) 1 ^ is the geometrical average of the three harmonic frequencies. Notice that in this limit the 
system, differently from traditional thermodynamic bodies, is characterized by a non uniform density varying on a 
macroscopic scale, much larger than the typical interatomic distance. The critical temperature takes the value 

k B T c = O.MhuhoN 1 / 3 (3) 



while the condensate fraction obeys the law 



for T <T C . Similarly to the case of the uniform gas one predicts a phase transition even in the absence of two-body 
interactions. This unique prediction of quantum statistical mechanics can be now tested experimentally in trapped 
Bose gases ffijl]. In fact, close to the critical temperature, these gases are so dilute that one can neglect in first 
approximation the role of the interaction. Estimate (^) has actually guided experimentalists to the right domain of 
temperatures which turns out to be of the order of fj,K (typically the trapping frequencies are of the order of a few nK 
while the number N of atoms varies from 10 4 to 10 7 ). Also the temperature dependence of the condensate fraction 
turns out to be in reasonable agreement with the theoretical prediction (||) as shown by fig 1. 

An interesting question is how two-body forces modify the ideal gas estimate (Q). Interaction effects can have two 
different origins. A first one is a genuine many-body effect exhibited by uniform systems where one works at fixed 
density. This effect is highly non trivial from the theoretical point of view but the resulting correction is expected 
to be small and difficult to observe in trapped gases. A second effect is the result of mean field forces which, in the 
case of positive scattering lengths, tend to push the gas towards the external region, thereby reducing the density 
of the gas. This causes a natural decrease of the critical temperature. The opposite happens if the interaction is 
attractive. The mean field effect can be easily calculated and the result takes the form || 

^ = -1.3-2-^/" (5) 



where a is the scattering length and aho = \/^/ rnu; ' 10 is the oscillator length. Typical values for 5T C /T C are a few 
percent in the available traps and the present accuracy of experimental data does not allow for a quantitative check 
of this prediction. By tuning the value of the scattering length to larger values it should be possible to obtain more 
sizable corrections to the value of the critical temperature. 

The above discussion has introduced a first important combination of the relevant parameters of the problem, given 
by N x / e a/ a} lo . This combination fixes the importance of interaction effects on the thermodynamic properties of the 
system. Its origin is easily understood if one calculates the ratio between the interaction energy Ei nt ~ gn where 
g = 4:Trh 2 a/m is the coupling constant fixed by the scattering length and n is the density of the thermal cloud evaluated 
in the center of the trap, and the thermal energy Et ~ ksT. A simple estimate for the thermal density is given 
by the classical value n ~ N(ksT '/muohof 1 ^ 2 ■ Using result (||) for the critical temperature, one finds that the ratio 
Eint I ' Et, evaluated at T ~ T c , is actually fixed by the combination 7V 1 / 6 a/a? lo . This should not be confused with the 
so called Thomas- Fermi combination Na/aho which instead characterizes the effects of two body- interactions on the 
zero temperature properties of a trapped condensate. The Thomas-Fermi combination is recovered by evaluating the 
ratio between the interaction energy Ei nt ~ gn, where n is now the central value of the condensate density, and the 
quantum oscillator energy Tuvq. The dependence on N given by the two combinations discussed above is extremely 
different and explains why the effects of the interaction at T ~ T c are small, while at T ~ they are so large that 
cannot be treated in a perturbative way. The role of the Thomas-Fermi parameter will be discussed systematically 
in the next sections. 

Before concluding this section it is useful to remind that the measurements of the temperature dependence of the 
condensate density available in trapped Bose gases are more reliable than the ones obtained in superfluid helium M 
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(see fig. 2). In helium the condensate fraction is not a natural observable and can be extracted only through elaborated 
analysis of neutron scattering data at high momentum transfer. Actually the main evidence for the phase transition 
in superfluid helium comes from the analysis of other thermodynamic quantities, like the viscosity and the specific 
heat. 



C. Order parameter and long range order 

A key quantity characterizing the long range order of a many-body system is the one-body density matrix M : 

n(r,r') = (* f (r)*(r')) (6) 

where ^(r) = ao</>o( r ) + J2i^o l ^'M r ) ^ s the field operator expressed in terms of the annihilation operators relative 
to a generic basis of single particle wave functions. We have here separated the contribution arising from the lowest 
single particle state <fio in order to emphasize the effect of BEC. 

Let us first consider an ideal gas at zero temperature. In this case all the particles occupy the same state $0 
determined by the solution of the Schrodinger equation for the single particle Hamiltonian — ?i 2 V 2 /2m + V ex t- The 
density matrix (||) then takes the separable form n^(r,r') = A</>Q(r)</>o(r') and remains different from zero for 
macroscopic distances | r — r' | of the order of the size of the sample (long range order). This is also called first order 
coherence and is a key consequence of Bose statistics. 

At finite temperatures one should include the thermal occupation of the other single particle states and eq. @ yields 

ra(r,r') = No^M?') + E M>*(r)<Mr') (7) 

where rii = (exp /3(ej — fi) — l) -1 are the thermal occupation numbers and a are the single-particle excitation energies of 
the excited states <pi . In the thermodynamic limit the ratio Nq /N tends to a constant value (condensate fraction) giving 
rise to long range order, while the second term in the r.h.s. of eq.(^) can be replaced with an integral and consequently 
vanishes at large distances. The importance of long range order is hence reduced at finite temperature since N is 
smaller than at T — 0. For a uniform ideal gas the decay of the thermal component is given by n(s) — > kTm/ [Airs) 
where s =| r — r' |. Above the critical temperature one has no more long range order (Nq/N — ► 0) and the one-body 
density matrix decays faster than 1/s at large distances. In the limit T T c one finds the gaussian behaviour 
n(s) = n(0)exp[-s 2 mk B T/2h 2 }. 

The effects due to Bose-Einstein condensation discussed above can be also derived by setting So = = \/~No. This 
is the well known Bogoliubov prescription which corresponds to ignoring the non commutativity between the particle 
operators a,Q and and treating them like c-numbers. In virtue of the Bogoliubov assumption the field operator 
takes the form 

*(r) = * (r) + 5><Mr) (8) 



i#0 



where 



*o(r) = VNoMr) (9) 



is a classical field, called the order parameter. The order parameter can be also regarded as the expectation value of 
the field operator 

*o(r) = (*(r)) (10) 

where the average is taken on a configuration with broken gauge symmetry. In terms of the order parameter the 
diagonal term of the density matrix can be written as 

n(r) = n (r) + n T (r) (11) 

where no(r) =| \l/o( r ) | 2 is the condensate density and riT(r) = ni I <M r ) 1 2 i s the contribution arising from the 

particles out of the condensate. In a dilute gas Ut can be safely identified with the density of the thermal component. 
In general one should not however ignore the fact that even at zero temperature the total density is different from the 
condensate density since the occupation numbers rii with i ^ of eq. (M) differ from zero, giving rise to the quantum 
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depletion of the condensate. An example where quantum depletion plays a crucial role is presented in fig. 3 which 
reports the density profile of a cluster of helium atoms calculated at T — using a correlated basis approach [||. 
These systems, differently from atomic gases, are strongly correlated. The figure clearly shows that the contribution 
of the condensate is only a small fraction (about 10%) and that in this case the measurement of the density profile 
would not yield any useful information on the condensate density no even at zero temperature. In the calculation 
of ref. [^) the condensate density has been evaluated by determining the natural orbits which diagonalize the 1-body 
density matrix (see eq.(0)) 

D. Equation for the order parameter 

The equation for the order parameter can be derived starting from the Heisenberg equation for the field operator 







ih—V(r,t) = [*(r,t),H] = - — + V ext (r,t) 



dt 



2yj2 



2m 



J &(r',t)V(r' -r)V(r',t)dr' $(r,t) . (12) 



This is an exact equation if one uses for V(s) the exact two-body interaction between particles and holds not only for 
dilute gases, but also for strongly interacting systems like superfluid 4 He. A closed equation for the order parameter 
is obtained if we replace the field operator ^ with the classical field Wo or i equivalently, with its expectation value 
(|Io|). This procedure is correct provided a series of assumptions are satisfied: 

i) The thermal depletion should be small so that the thermal fluctuations of the field operator can be ignored. This 
implies temperatures much smaller than the critical value. 

ii) The system should be so dilute that the quantum fluctuations of the field operator can be ignored. This requires 
that the gas parameter na 3 be much smaller than unity or, equivalently, that the average distance between particles 
be much larger than the scattering length. 

hi) The microscopic interaction V(s) should be replaced by the pseudopotential gAirTi 2 {a/m)5(s) fixed by the s-wave 
scattering length a. This implies that only low energy features of the problem can be investigated with this approach. 
Equivalently, only variations of the order parameter over distances larger than the range of the interaction can be 
explored. 

iv) The number of atoms N should be large enough in order to justify the Bogoliubov prescription. 
Under the hypothesis i)-iv) one obtains the most famous Gross-Pitaevskii equation || 

/ ft 2 v 2 \ 

ih— * (r, t) = I + V ext (v, t)+g\ * (r, t) | 2 J tf (r, t) , (13) 

for the classical field This equation provides the basic tool to explore the static and dynamic features of non 
uniform Bose gases at low temperature. The GP equation shares interesting analogies with the Maxwell equations 
of classical electromagnetism. In both approaches one provides a description of a many-body problem (atoms anf 
photons respectively) in terms of a classical field (order parameter and electromagnetic field ) . From this point of view 
a gas of Bose-Einstein condensed atoms can be regarded as a classical matter wave. With respect to the equations 
of electromagnetism the GP equation however exhibits important differences: first it contains an important non linear 
term due to interactions (this difference is not crucial since non lineariteies are exhibited also by electromagnetic 
phenomena in dispersive media); a second important difference is the role played by the Planck constant Ti. Differently 
from the Maxwell equations the GP equation (|l^) in fact depends explicitly on this fundamental constant of quantum 
mechanics. This feature is directly connected with the different momentum-energy relationship exhibited in the two 
cases. For massive particles the relation is quadratic (E = p 2 /2m), while for photons it is linear (E — cp). When we 
employ the ondulatory description through the de Broglie prescription E = Tim and p — hk, the relationship bewteen 
the wavevector and the frequency is still independent on the Planck constant in the case case of photons (u> — ck), 
while it exhibits an explict dependence on Ti in the case of atoms (w = fik 2 /2m). This explains, for example, why 
interference phenomena with classical matter waves depend explicitly on the value of the Planck constant. 



E. Ground state configuration 



If we look for stationary solutions of the form 
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*o (r, t) = * (r) exp(-v*t/ft) (14) 
the Gross-Pitaevskii equation ([l3]) becomes 

?i 2 V 2 



. !(; +F ext (r)- M + g|* (r)|^*o(r) = 0. (15) 

This equation can be also derived using the standard variational procedure 

5(E-fiN) = (16) 

where 

E = J 1 V *° |2 I *o I 2 +| I *o I 4 ) <*r - £ fcm + £ eKt + E mt (17) 

is the energy functional of the system. In this equation we have naturally identify the kinetic energy (Eki n )> the 
external potential (E ext ) and the mean field interaction (Ei nt ) energies. Starting from the GP equation one can easily 
derive the useful relationship 

N[i = E km + E ext + 2E mt (18) 
for the chemical potential, and the virial identity 

2Ekin — 2E ext + 3E int = (19) 
characterizing the equilibrium configuration. It is also useful to recall the so called release energy 

-^release — -^kin ^int (^0) 

which coincides with the energy of the system after switching off the trap and is hence given by the sum of the kinetic 
and interaction energies. 

It is important to remind here that the quantity Eki n entering the above equations is the kinetic energy of the 
condensate (mean field kinetic energy) and should not be confused with the full kinetic energy of the system. The 
latter takes in fact contribution also from the particles out of the condensate and, using the notation (g) for the 
1-body density matrix, can be written in the form 

7 total _ f | I v-.t. 2 \T 



E k?n = I dr— | V* (r) \ 2 +^n t | V&(r) \ 2 ] . (21) 
"' 

While for dilute gases the contribution n * °^ ^ ne 11011 condensate component to the total number of atoms is small 
at T — and one can safely identify N with AT , the corresponding contribution to the kinetic energy (second term of 
eq.(|2l])) may be much larger than the condensate kinetic energy. For similar reasons the mean field interaction energy 
Ei nt should not be confused with the full interaction energy whose evaluation requires the explicit knowledge of short 
range correlations. The distinction between the mean field and total energies becomes particularly clear and stark in 
the case of a uniform gas interacting with hard sphere potentials. In this case the total kinetic energy coincides with the 
energy of the system, the interaction energy being exactly zero due to the special choice of the potential. Conversely 
the kinetic energy of the condensate is zero because of the homogeneity of the sample. Notice that the mean field 
scheme fails to evaluate the full kinetic energy ( pl| ) . This can be seen in the case of a uniform gas where the occupation 
number of the single particle states is given by the Bogoliubov expression n p = (p 2 /2m + mc 2 )/2e(p) — 1/2. The 
quantity Xp^o n v 1S evaluable and gives the quantum depletion of the condensate. Viceversa the sum ^2i^ {p 2 /2m)n p 
exhibits an ultraviolet divergency, showing that the evaluation of the kinetic energy requires the proper knowledge of 
the momentum distribution at momenta of the order of the inverse of the scattering length and cannot be achieved 
within the mean field scheme. 

It is finally worth pointing out the key role played by the chemical potential in the formalism of the GP equation. 
As one can see from ( |T4| ) the time dependence of the order parameter is fixed by fi and not by the energy. This 
is a consequence of the fact that the order parameter is not a wave function and that the GP equation is not a 
Schrodingcr equation in the usual sense of quantum mechanics. From the point of view of many-body theory the order 
parameter corresponds to the matrix element of the field operator between two many-body wave functions containing, 
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respectively, N and N + 1 particles. This implies that its time dependence is fixed by the factor e -*( E (-' v + 1 )- E W)* 
and hence by the chemical potential // = dE/dN rather than by the energy E. 

If the Thomas Fermi parameter Na/aho is very large the GP equation ([l3|) can be solved analytically. In fact in 
this case the density profile, as a consequence of the repulsive effect of the interactions, becomes so smooth that the 
kinetic energy term in dUf ) can be ignored and the density takes the shape of an inverted parabola (we consider here 
spherical trapping) 

n(r) = i(/,-F eKt ) = ^(l--^) (22) 

where the radius R of the condensate is given by the formula 

/ \ 1/5 

R = a ho 15N— (23) 
V a ho J 

and /i = mui 2 R 2 /2. For an axially deformed trap the density profile is simply generalized to n(y) — (/// ' g){\ — r\_j 'R\ — 
z 2 /Z 2 ) where r 2 ^ — x 2 + y 2 and the radial and axial radii are related by the equation muj 2 ± R 2 L = m,io 2 z Z 2 — 2/i. In 
many available configurations the size of the atomic cloud can become much larger than the oscillator length and 
reach almost macroscopic values. 



F. Visibility of the condensate 



In the introduction we have emphasized the important fact that, due to the harmonic trapping, the condensate 
density can be separated from the thermal component, a crucial condition for the experimental visibility of the 
condensate. Due to the importance of this fact it is worth discussing in a deeper way the consequences of two-body 
interactions. In fact the argument usually invoqued is based on the ideal gas model where the width of the condensate 
is fixed by the oscillator length aho — y/tt/muho- On the other hand the width of the thermal cloud is of the order of 

, rp \ 1/2 



Rt „ - (24) 

so that, for "macroscopic" temperatures fcsT> Tiujho, one has Rt ^> aho with the consequent net separability of the 
two components. In the previous section we have however shown that interactions are responsible for a huge increase 
of the size of the condensate whose value, in the Thomas- Fermi limit, is given by eq.(p3|). One then conclude that 
interacations will partially reduce the visibility of the condensate by making its size comparable to the one of the 
thermal cloud. Interactions, by producing a smoother density profile, make the system closer to a uniform gas and 
consequently reduce the visibility of BEC in coordinate space. In fig 4 we show a typical in situ measurement of 
the density profile taken at different temperatures JlO[ . A similar argument holds if one explores the behaviour of 
the release energy E re i ease . For an ideal gas the release energy is fixed by the oscillator energy Htoo, which is much 
smaller than the release energy of the thermal cloud, fixed by ksT. This provides a net separation between the two 
components which expand at different velocities. The difference is sizably reduced if one takes into account interactions 
which increase significantly the value of E re i ease making it comparable to the thermal energy kT. In conclusion the 
separability of the condensate from the thermal component is significantly reduced by two-body interactions and 
requires a careful bimodal fit to the measured profiles (both in in situ and after expansion) . 

While the visibility of the condensate in coordinate space is reduced by two-body (repulsive) interactions, it is 
conversely enhanced in momentum space. In fact, according to the Heisenberg relationship, the width of the condensate 
in momentum space behaves like h/R and hence becomes smaller and smaller as R increases, approaching the typical 
(5 distribution characterizing BEC in uniform systems. While in the absence of interactions the visibility of BEC 
in coordinate and in momentum space are perfectly equivalent, due to the symmetric role played by the spatial and 
momentum coordinates of the harmonic Hamiltonian, in the Thomas- Fermi limit this symmetry is lost. The possibility 
of measuring the momentum distribution of trapped atomic clouds, via inelastic photon scattering experiments |Tl] ] 
opens new possibilities to separate at a deeper level the condensate from the thermal component. In principle such 
measurements might provide an identification of the quantum depleted part too. In practice the effect of quantum 
depletion is too small to be observed in the presently available condensates. 



G 



G. Role of the phase: interference with BEC 



The phase of the order parameter plays a crucial role in characterizing the interference phenomena exhibited by a 
Bose-Einstein condensed gas. The simplest example is the interference produced by two initially separated condensates 
which expand and overlap. The corresponding experiment was first carried out at Mit 

Let us first consider a single expanding condensate. At large times the phase of the order parameter takes a 
quadratic dependence, yielding the asymptotic form 

*o(r,t)=|*o(r,t) \ e ^/2M (25) 

for the order parameter. It is interesting to notice that the velocity field 

v = -VS (26) 

TO 

associated with the order parameter (^) coincides with the classical law v = r/t. For large times the velocity then 
exhibits a fully classical behaviour. However the phase S keeps its quantum nature fixed by the Planck constant. 

If we now consider two expanding condensates separated by a distance d and we assume that for large times the 
order parameter is given by a linear combination of the form 

*o =| *i | e "™ (r+d/2)2/(2?it) + e ia | * 2 | e im{r ~ A/2)2 /(2M) (27) 

we find that the density | | 2 acquires a modulation with typical fringes orthogonal to the vector d and separated 
by the distance 

< 28 > 

This distance depends explictly on the Planck constant and with typical choices of the parameters d and t is a visibile 
length of the order of 10 — 20(iK. The above estimate for the interference effects ignores the role of the interaction 
between the two clouds which are responsible, during the expansion, for deviations from the law (|2?]). The basic 
physics is however accounted by this simple model which provides a useful qualitative estimate of interference effects. 
It points out, in particular, a conceptually important question connected with the role of the relative phase of the two 
condensates. Suppose that the two condensates are built in an independent way before measurement. Why should 
the order parameter have the form ( p7j ) with a well defined value of the relative phase a giving rise to interference? 
The question addresses the important problem of the quantum measurement and of the consequent reduction of the 
wave packet. The standard theoretical point view, based on the traditional rules of quantum mechanics, predicts that 
after measurement the wave function of the whole system will be of the form ( p7j ) with a relative phase which cannot 
predicted in advance, unless the condensates were already in a coherent configuration before the expansion. 



H. Role of the phase: irrotationality of the superfluid flow 

A crucial consequence of the existence of an order parameter of the form (1), is the irrotationality of the flow. This 
feature is expressed by eq.(^6|) from which one deduces that the phase plays the role of a velocity potential. Actually 
starting from the Gross-Pitaevskii equation for the order parameter it is possible to derive the equations of motion in 
a form which resembles the equations of irrotational hydrodynamics. These are the basic equations for the describing 
the motion of a superfluid. The derivation is straightforward in a dilute gas. One starts from expression (1) for the 
order parameter and derives, using the GP equation a closed set of equations for the density and for the phase 
or, equivalently, for the velocity field. The equations take the form ]l3[ ] 

^-+div(vn) = 0, (29) 



9 f 1 Ti \ 

h 7^ S + o mv2 + v ** + -9" ~ o 7= y2 V™ ) = 0- (30) 

ot \2 2myn J 

The first equation is the usual equation of continuity, while the second one can be regarded as an equation for 
the velocity potential. Notice that in these equations, which are exactly equivalent to the original Gross-Pitaevskii 



7 



equation, the Planck constant enters through the so called quantum kinetic pressure term proportional to ~S7 2 y/n. 
This term should be compared with the interaction gn given by two-body forces. If the role of interactions is very 
important the quantum pressure term can be neglected. This happens in the study of the equilibrium profile in the 
Thomas-Fermi limit. The quantum pressure term is negligible also in the time dependent problem, provided we study 
phenomena where the density varies in space on lengths scales L larger than the so called healing length 



^ y 2mgn ^ ^ 

If L ^> £ it is immediate to see that the quantum pressure in (|3(]) is negligible and that the same equation reduces to 
the simpler Euler-like form 

m + V Qtov 2 + V ext + = (32) 
where n — gn is the chemical potential evaluated for a uniform body at the corresponding density. The equation of 



continuity and the Euler equation (32) are called the hydrodynamic equations of superfluids. They are characterized 
by two important features: the irrotationality of the motion and the crucial role of interactions which suppress 
the effects of the quantum pressure. It is worth noticing that the " hydrodyinamic" form of these equations is not the 
result of collisional processes as happens in classical gases, but is the consequence of superfluidity. 

It is also useful to recall that the validity of the hydrodynamic equations is limited to the low T regime (at finite 
temperature one should use the formalism of two- fluid hydrodynamics). In the limit of small temperatures these 
equations can be used to study also dense superfluids, like 4 He as well as Fermi superfluids. Of course in this case 
one should use the corresponding expression for the chemical potential. It is finally worth mentioning that in the 
hydrodynamic limit the Planck constant, which entered equation ( |30| ) for the velocity field through the quantum 
pressure term, has disappeared from the equations of motion. 

Starting from the hydrodynamic equations of superfluids one can explore different interesting problems in trapped 
BEC gases. These include the propagation of sound and the excitation of collective modes. 



I. Sound and collective excitations 

In the limit of small perturbations (linear limit) the hydrodynamic equations ( p9[ ) and (|3^) can be reduced to the 
typical form 

^5n = V • [c 2 (r)VH (33) 

characterizing the propagation of sound waves in nonuniform media. Here c(r) is a r-dependent sound velocity fixed 
by the relationship mc 2 (r) = u — V^^r). Sound can propagate in trapped Bose gases if the wavelength A of the wave 
is smaller than the size of the system. In the Thomas-Fermi limit this condition can be well satisfied especially in the 
elongated direction of the condensate. In the experiments carried out at MIT |l4|] the axial size Z is a few hundred 
micrometers while £ is a fraction of micrometer, so there is wide space for the propagation of sound waves in such 
samples. It is also interesting to consider the case where the wave length A is smaller than the axial size (A < Z), but 
larger than the radial size (A > i?j_). In this case the nature of the wave is characterized by typical 1-dimensional 
features p|. 

If the wavelength becomes comparable to the size of the condensate the solutions of eq. (]33|) cannot be described in 
terms of a localized propagation of sound, but should be determined globally. These solutions correspond to oscillations 
of the whole system and the discretization of the eigenvalues cannot be ignored in this case. For spherical trapping 
one can obtain analytic solutions of ( |33"| ) in the form 5n — R nr iYim where I = 0,1,2.. is the angular momentum 
(in units of K) carried by the excitation, n r is the number of the radial nodes and R n e is the radial function to be 
determined by solving the equations of hydrodynamics. The corresponding eigenfrequencies obey the dispersion law 




^2n 2 r + 2n r i + 3n r + i (34) 

which should be compared with the prediction w^ (2n r + I) of the ideal gas model. The role of interactions is 
particularly evident for the surface excitations (n r — 0). In this case the hydrodynamic theory gives Vlujho to be 
compared with value £ojq of the ideal gas prediction. 
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Result (gj) reveals that, in the Thomas- Fermi limit Na/a^o 3> 1, the dispersion relation of the normal modes of 
the condensate has changed significantly from the noninteracting behavior, as a consequence of two-body interactions. 
However it might appear surprising that in this limit the dispersion does not depend any more on the value of 
the interaction parameter a. This differs from the uniform case where the dispersion law, in the corresponding 
phonon regime, is given by lo = cq and depends explicitly on the interaction through the velocity of sound. The 
behavior exibited in the harmonic trap is well understood if one notes that the discretized values of q are fixed 
by the boundary and vary as 1/R where R is the size of the system. While in the box this size is fixed, in the 
case of harmonic confinement it increases with N due to the repulsive effect of two-body interactions (see eq.(p3|): 
R ~ (A/ r a/aho) 2 ^ 5 (wwho) -1 ^ 2 - On the other hand the value of the sound velocity, calculated at the center of the trap, 
is given by c = (Naj ' a\ LO ) 2 ^{w\ lo /'m) 1 / 2 and also increases with N. One finally finds that in the product cq both the 
interaction parameter and the number of atoms cancel out, so that the frequency turns out to be proportional to the 
bare oscillator frequency Who- A similar argument holds also for the surface excitations. In fact in the presence of an 
external force the dispersion relation of the surface modes obeys the classical law J 1 = Fq/m where F — vhlo\ R is 
the force evaluated at the surface of the system. Since the product hqR gives the angular momentum carried by the 
surface wave, one immediately recovers the dispersion law Viuiho- 

Analytic solutions of the HD equation (33) are available also in the case of axi-symmetric potentials of the form 
Vext = ( m /^)( LU ± r '± + w z z2 )- I n this case the third component of angular momentum is still a good quantum number. 
For very elongated traps (uj z <C ui±) the frequency of the lowest solution of even parity takes the value Jlj| sfhj2u) z 
in excellent agreement with the experimental results of fll6| . 

The high precision of frequency measurements is not only providing us with a powerful tool to check the predictions 
of the hydrodynamic theory of superfluids, but can be also used to explore finer effects, like the temperature dependence 
of the collective frequencies and, possibly, beyond mean field corrections. 



J. Moment of inertia 



An important consequence of the hydrodynamic theory of superfluids is that the response 

9 = ltmn-,0 < o > ( 35 ) 

to a rotating field of the form H rot = —£IL Z , where L z — ~i^^2 k Yk x is the third component of the angular 
momentum operator, is given by the irrotational value 

e=f ^4^4^ ) @ ri 9 (36) 

of the moment of inertia where & r i g is the classical rigid value. This result can be easily derived by rewriting the 
equations of hydrodynamics in the frame rotating with the angular velocity fl. Eq.(|36|) can be also written in the 
form = e 2 ©i rr where e = (uj 2 — oj 2 )/uj 2 + uj 2 ) is the deformation of the trap in the plane of rotation and we have 
used the Thomas- Fermi results for < x 2 > and < y 2 >. 

The quenching of the moment of inertia with respect to the rigid value was confirmed in a series of experiments 
of the 60 's caried out on superfluid helium, providing an independent measurement of the superfluid density p7| . In 
this context it is worth mentioning that the superfluid density of superfluid helium exhibits a very different behaviour 
with respect to the condensate density, especially at low temperature. In particular at T = the superfluid density 
coincides with the total density while the condensate density is only a small fraction (~ 10%). 

A challenging question is how to measure the moment of inertia of a trapped gas where the direct measurement of 
the angular momentum L z is not feasible due to the dilutcness of the sample. An interesting possibility is provided by 
the fact that, if the deformation of the trap is different from zero, the quadrupole and rotational degrees of freedom 
are coupled each other. This is well understood by considering the exact commutation relation 

[H,L z ]=im(u J 2 x -u 2 y )Q (37) 

which explicitly points out the link between the angular momentum operator L z and the quadrupole operator Q = 
~Ylii x iVi- Since the quadrupole variable can be easily excited and imaged one expects to obtain information on the 
rotational properties of the system by investigating the quadrupole modes excited by the operator Q. An important 
case is the so called scissors mode. This mode corresponds to the rotation of a deformed condensate whose inclination 
angle 9 oscillates in time. The oscillation can be induced by a sudden rotation of the confining trap with respect to 
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the initial equilibrium configuration. It is important to point out that the frequency of this oscillation does not vanish 
with e as one would expect for a classical system. In fact in a superfluid both the restoring force parameter and the 
mass parameter (moment of inertia) behave like e 2 so that the frequency keeps a constant value when e — > 0. The 
hydrodynamic theory of superfluids predicts the value p8| 

uj hd = y^ 2 + w 2 (38) 

in excellent agreement with the experimental results recently carried out at Oxford Jl9[ ] (see fig. 5). Conversely in 
a classical gas one predicts the value ui =| uj x ± ui y |. Also these values have been tested experimentally with high 
accuracy by exciting the rotation of the gas at high temperatures, well above T c . 



K. Quantized vortices 

Quantized vortices are one of the most spectacular manifestations of superfluidity (we will not discuss here other 
important manifestations of superfluidity, like the reduction of viscosity recently explored at MIT |2(J ) . The existence 
of quantized vortices is the combined consequence of the behaviour of the phase of the order parameter, which 
fixes the irrotationality of the velocity field, and of the non-linearity of the equations of motion which is a crucial 
the consequence of two-body interactions. In the Gross-Pitaevskii theory a quantized vortex can be regarded as a 
stationary solution of eq. ( |l3|) of the form 

*o(r,t) = e- i ' lt e^^(r) (39) 
where cj> is the azimuthal angle and *Sf v is a real function obeying the equation 
ft 2 V 2 



2m +^+2^ + U &)+9*i(r ± ,z) 



$v(r_L,z)= fx^ v (r x ,z) . (40) 

The velocity field associated with the order parameter ( |39"| ) takes the form v = (fi,/m)V0 = Ct x r/r 2 where $7 is the 
unit vector along the 2-th direction. It satisfies the irrotationality constraint everywhere except along the vortical line 
and gives rise to a total angular momentum given by 

L z = Nh . (41) 

The " centrifugal" term in 1 jr\ of eq. (f40| ) originates from the peculiar behaviour of the velocity field of the vortical 
configuration. This term is responsible for the vanishing of the condensate density | ^ v | 2 along the z-th axis and 
produces the so called vortex core, whose size is fixed by the healing length. In fig. 6 we plot a typical vortical profile 
obtained by solving the Gross-Pitaevskii equation (f4(i|). 

Quantized vortices are excited states of the system but they should not be confused with the elementary excitations 
discussed in sect. I . From the many-body point of view the vortex correspond to a state where all the N particles of 
the sample occupy the new solution of the GP equation. The excitation energy and the angular momentum associated 
with the vortex are consequently " macroscopic" quantities while the excitation energy and the angular momentum of 
elementary excitations are of the order, respectively, of the trapping frequencies and of a few units of fi. 

The excitation energy E v of the vortex, given by the difference between the vortex and the ground state energies, 
can be calculated starting from the solution of the GP equation ( |40| ) a nd using expression (17) for the energy. For 
large N samples the calculation can be done analytically. One finds plf 

4tt h 2 , 0.67i? ± 

E v = —n Z— In —± (42) 

3 m £ 

where uq is the value of the density in the center of the trap. It is worth noticing that the factor in front of the 
logarithm is proportional to the so called column density J dzn(x,y, z) — (A/3)Zno calculated along the symmetry 
axis (x = y = 0). As we will see later this identification is useful for the calculation of the the energy of a vortex 



line displaced from the symmetry axis. Result 
generate an energetically stable vortex. In fact, ca 



allows for an estimate of the critical angular velocity needed to 
Tculating the energy H — QL Z of the sample in the frame rotating 
with angular velocity f2, one easily finds that the vortex is the lowest energy configuration if the angular velocity f2 
exceeds the critical value 
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E v 5 h 0.67R ± 
^ = LT = 2^i ln ^— ' (43) 

In deriving result ( ^3|) we have used the relationship N — (8/15)7rnoR 2 L Z . Typical values of Q cr correspond to fractions 
of f2j_, the exact value depending on the geometr y o f the trap and on the value of the Thomas-Fermi parameter. In 
the recent experiments carried out at ENS-Paris [£2| with an enolangated trap (uj z /lj±_ ~ 1/10) the calculated value 
of Q C r turns out to be about oj_l/3. In this experiment the vortex is created by rotating a slightly asymmetric trap 
at different values f2, in analogy with the rotating bucket experiment on superfluid helium. This procedure differes 
from the one followed at JILA ]2j| where the vortex was instead created with full optical methods. The value of the 
critical angular velocity measured at ENS is a factor 2 higher than the theoretical estimate (fl3]). The discrepancy is 
due to the fact that, even if the vortex is energetically favourable, its nucleation is forbiddeen by the occurrence of 
a barrier which disappears only at higher angular velocities. The nucleation of the vortex as well as the study of its 
stability and life time is at present a topic of intense experimental and theoretical investigation. 

An important question in the study of vortices is not only how to produce but also how to detect them. The size 
of the vortex core is in fact too small (fractions of microns) to be observed in situ. It can be however observed by 
imaging the atomic cloud after expansion. An alternative method is given by the measurement of angular momentum, 
a quantity which, according to eq. ( ^jj ) , characterizes in a peculiar way the quantum nature of the vortex. The presence 
of angular momentum produces a splitting of the quadrupole frequencies tu± relative to the excitations with £ z = ±2 
( £ z is the third component of the agular momentum of the elementary excitation) . The value of the splitting and its 
relation with the angular momentum L z of the vortex can be derived using a sum rule approach. The result, in the 
Thomas- Fermi limit, is given by the formula p3| 

lo+ - w_ = 2 <L % > . (44) 

The splitting of the quadrupole frequencies lu± produces a mechanism of precession that can be detected experimen- 
tally. To this purpose one generates a small quadrupole deformation in the plane orthogonal to the vortex axis, by 
adding a field of the form x 2 — y 2 . The precession of the quadrupole deformation is directly related to the splitting 
©. In fact, by requiring that in the frame rotating with angular velocity f2 prec the energies lu± of the two modes 
are degenerate, one finds the result fl pre c = { UJ + — w -)/4- Using the above procedure the authors of ref p4j have 
observed a jump in the angular momentum of the condensate above a critical value of the stirring frequency (see fig. 
7). This approach has permitted to test the quantization of the angular momentum of a single vortex line. 

Gyroscopic effects of different type have been observed in [BBj where the precession of a vortex line circulating 
around the symmetry axis of the condensate was investigated. The angular velocity characterizing this precession 
can be estimated by evaluating the energy and the angular momentum of a vortex line as a function of its distance d 
from the symmetry axis. This calculation is easy if one assumes that the vortex line is straight. The d dependence 
of the energy is determined, within logarithmic accuracy, by the c?-dependence of the column density J drn(x,y, z) 
characterizing the prefactor of eq.(^). This yields [?] 

/ d 2 \ 3/2 

E v = Nhn c 1^1 - w ) (45) 

where fl c is the critical angular velocity (|43j). Using Stokes' theorem and working in the Thomas- Fermi limit it is 
possible to calculate also the angular momentum for which one finds the result 

/ d 2 \ 5/2 

L z = Nh(l- w ) . (46) 

The angular velocity characterizing the precession is given by fl prec = dE v /dL z . For small values of the vortex dis- 
placement (d <C R±) one obtains the result 26 0, prec = (3/5)f2 cr which agrees reasonably well with the experimental 
results of JILA §7\. 



L. Conclusions and acknowledgements 



The main purpose of these notes was to give an introduction to some key features which have been the object of 
recent experimental and theoretical work in the physics of trapped Bose-Einstein condensates. Special emphasis was 
given to the interplay between the phase of the order parameter and two-body interactions which characterize the 
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non linearity of the Gross-Pitaevskii equation. This interplay is at the basis of the superfluid behaviour exhibited by 
these systems and is responsible, in particular, for the existence of quantized vortices. 

It is a pleasure to thank the fruitful collaboration of the BEC Trento team and the warm hospitality of the French 
friends at the Institute d'Etudes Sientifiques de Cargese. 
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FIGURE CAPTION: 

Fig.l. Condensate fraction of a trapped Bose gas as a function of T/T c . Circles are the experimental results of [||, 
while the dashed line is eq.([|). 

Fig. 2. Condensate fraction as a function of temperature in superfluid 4 He (from |(|). 

Fig. 3. Total density and condensate density profiles of a helium droplet of 70 atoms (from Q). 

Fig. 4. Axial profiles of a cloud of sodium atoms at different temperatures (from JhJ). 

Fig. 5. Scissors oscillation observed in JljJ. The full line corresponds to a sinusoidal fit with u/2tt = 265 ± 0.8 Hz to 
be compared with the theoretical prediction 265 Hz of eq.|38|. 
Fig. 6. Density profile of a quantized vortex (from 0). 

Fig. 7. Angular momentum per particle of a trapped condensate as a function of the angular velocity of the trap (from 
©)■ 
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